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ABSTRACT 



Heyer , and Michael G. Burton 



The Galactic supernova remnant (SNR) IC 443 is one of the most studied core-collapse SNRs for 
its interaction with molecular clouds. However, the ambient molecular clouds with which IC 443 
is interacting have not been thoroughly studied and remain poorly understood. Using Five College 
Radio Astronomy Observatory 14m telescope, we obtained fully sampled maps of ~ 1° x 1° region 
toward IC 443 in the 12 CO J = 1 — and HCO + J = 1 — lines. In addition to the previously known 
molecular clouds in the velocity range ^lsr = — 6 to —1 krns" 1 (—3 kms -1 clouds), our observations 
reveal two new ambient molecular cloud components: small (~ 1') bright clouds in ^lsr = — 8 to 
—3 kms - (SCs), and diffuse clouds in ^lsr = +3 to +10 kms - (+5 kms - clouds). Our data 
also reveal the detailed kinematics of the shocked molecular gas in IC 443, however the focus of this 
paper is the physical relationship between the shocked clumps and the ambient cloud components. 
We find strong evidence that the SCs are associated with the shocked clumps. This is supported by 
the positional coincidence of the SCs with shocked clumps and other tracers of shocks. Furthermore, 
the kinematic features of some shocked clumps suggest that these are the ablated material from the 
SCs upon the impact of the SNR shock. The SCs are interpreted as dense cores of parental molecular 
clouds that survived the destruction by the pre-supernova evolution of the progenitor star or its nearby 
stars. We propose that the expanding SNR shock is now impacting some of the remaining cores and 
the gas is being ablated and accelerated producing the shocked molecular gas. The morphology of 
the +5 kms -1 clouds suggests an association with IC 443. On the other hand, the —3 kms -1 clouds 
show no evidence for interaction. 

Subject headings: supernova remnants - ISM: individual (IC443) 
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1. INTRODUCTION 

Massive, newborn stars subject the surrounding parent 
molecular cloud to disruptive ultra violet radiation and 
stellar winds (e.g., iChevalierl Il999h . For O-type stars, 
a large area around them is expected to be cleared of 
molecular clouds, thus when the stars explode as super- 
novae (SNe), their remnants will spend a significant time 
expanding inside the hot bubble created by the progeni- 
tors. On the other hand, the destruction of natal molec- 
ular clouds would not be as effective for early B stars, 
and the supernova remnants (SNRs) of these less mas- 
sive stars could be able to interact with the molecular 
clouds that survived the destruction. Therefore, the en- 
vironment of core-collapse SNRs not only affects the evo- 
lution of SNRs, but also provides information about their 
progenitor star. 

Among the known 274 Galactic SNRs (|GreeiJl2QQ9h . a 
few tens of them show direct or indirect evi dence of in- 
terac tion with ambient molecular clouds (see Uiang et all 
1201(1 ). IC 443 (G189. 1+3.0) is one of the first and most 
studied SNRs for its interaction with molecular clouds. 
The overall morphology of IC 443 is depicted in Fig.[TJa). 
IC 443 has two shells with different radii (Shells A and 
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B, nomenclature is adopted from lBraun fc Stroml ll986). 
that are clearly seen in optical and radio continuum im- 
ages. This peculiar morphology has led to a suggestion 
that the r emnant in i tially evolved in an inhomogeneous 
medium. iLee et al.1 (|2008l LEE08 hereafter) suggested 
that Shell B could be a blowout part of the remnant into 
a rarefied medium. There is another larger and fainter 
sh ell (Shell C) that partia lly overlaps with Shells A and 
B. iBraun fc Stroml (|1986l ) proposed that Shell C is also 
physicall y associated with She l ls A a nd B, however, the 
work of lAsaoka fc Asc henbach (jl994j ) suggests that Shell 
C is more likely another SNR that is only positionally co- 
incident in the sky . 

IC 443 is apparently interacting with both atomic 
and molecular gas (LEE08, and references therein). 
The SNR shock in the northeastern part mostly shows 
the atomic lines expected in postshock recombining 
gas (Fes en fc Kirshnerl Qj)80), and the shock is prop- 
agating into an atomic medium. On the other 
hand, the southern part of the remnant shows shock- 
excited broad molecular lines, an d the SNR sho ck 
is encounter i ng m o lecular clouds ([ Denoyer 11979 blljl; 
iBurton e t al. 1988; Ziurvs et al. 1989; Dick man et al.l 
119921 : Ivan Dishoeck et al.l I1993D . The Two Micron Al l 
Sky Survey JHK composite image of iRho et al.1 (|200lh 
well demonstrates the different nature of the ambient 
medium. 

Shocked molecular gas associated with IC 443 was 
first di scovered in OH and 12 CO(J=l— 0) by iDenoyerl 
(|1979 b.a). Shock acceleration is directly indicated by the 
large velocity width of the observed lines (> 30 kms -1 ). 
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Figure 1. (a) Schematic showing the overall morphology of IC 4 43 adopted from LEE08, with the nomenclatures for shells overlaid. 
The asterisk indicates the position of the pulsar (Olbert et al. 2001). (b) 21cm radio continuum image of IC 443 (LEE08) with the area 
covered by the FCRAO observations: thin solid polygon for 12 CO (J=l— 0) and thick dotted polygon for HCO + and HCN. The contours 
are distribution of shocked 12 CO (see Fig.|4]for details). The shock clumps (B-H) identified by Denover (1979a) and Dickman et al. (1992) 
are indicated. 



Broad wings of various molecular lines have been re- 
ported (e.g. JZiurvs et al .11 19891: Ivan Dis hoeck et al. 1993) 
including H9O (jSnell et al.l l2Q05h . The overall distribu- 
tio n of shocked molecul ar gas arou nd IC 443 was mapped 
bv iBurton etHI (|1988ft in H 2 , bv lDickman et al.l (|1992l) 



in 12 CO ( J=l-0) and HCO+ ( J=l-0), and bv lXu et al.l 
(|20lTh in 12 CO (J=2-l, J=3-2). Shocked molecu- 
lar emission is bright along the southern boundary of 
Shell A, but is also found toward the northeast and near 
the center of Shell A (Fig. [TJb)). No shocked molecular 
gas outside Shell A has been reported. The overall distri- 
bution of shocked mol ecular gas has often been described 
as an expanding ring ( Burton et al.|[l98 8; Dick man et al.l 
119921 : Ivan Dishoeck et al.lll993ft . 

On the other hand, the identification and the nature 
of ambient clouds that are physically interacting with 
IC 443 have not been comprehensively studied. Ambient 
molecular clouds tow ard IC 443 were first reported by 
iCornett et al.l (j!977ft . They found a geometrically thin 
(~ 3 pc), sheet-like molecular clouds centered around 
f lsr ~ —3 km s" 1 and propo sed its interaction with 
IC 443. Dickman et "all ()1992h published more refined 
maps of the ambient molecular gas. IC 443 is located 
toward the Galactic anticenter direction where the ra- 
dial velocity due to the Galactic rotation is degenerate 
to f lsr °° Okms -1 . The clouds at ^lsr ~ — 3 kms -1 
partially overlap with IC 443 in the sky and they are of- 
ten assumed to be the clouds that are interacting with 
IC 443 (e.g. JTroja et al.ll2006ft . although no clear indi- 
cation of their physical association has been found. 

Identifying the ambient clouds that are physically asso- 
ciated with SNRs is important as it provides fundamental 
information about the environment of the SNRs. This 
is particularly important for core-collapse SNRs whose 
environment is significantly affected by stellar feedback. 
Also, the study of ambient clouds is important in under- 



standing the hadronic nature of associated 7-ray emis- 
sion. 7-ray emission can be emitted when cosmic rays 
that are accelerated at SNR shocks encounter nearby 
dense molecular clouds. IC 443 has long been suspected 
being a site of cosmic ray accele ration. It harbors a n 
EGRET source 3EG J0617+2238 (jHartman et al.l[l999h . 
which might be associate d with the hard X- ray source 
in the eastern boundary (jBvkov et all 12008ft . A TeV- 
source in the weste rn boundary is detected by MAGIC 
([Albert et al.ll2007ft and VERITAS (jAcciari et al.ll2009h . 
More recently, GeV e mission from IC 443 is re- 
ported by Fermi L AT (jAbdo et al.l 12010ft and AGILE 
(jTavani et aTll2010ft . 

In this paper, we report mapping observations of 
molecular lines toward the Galactic SNR IC 443. 
A square-degree field toward IC 443 is mapped in 
12 CO( J=l-0) and HCO+( J=l-0) with spatial resolution 
of ~ 50". Although IC 443 has been observed in various 
molecular lines, observations have been focused on the 
broad-line shocked molecular clumps. Our data present 
a global view of IC 443 and its interaction with molecular 
clouds. The 12 CO observations we report here were con- 
ducted using the Five College Radio Astronomy Obser- 
vatory ( FCRAO) 14m radio te le scope, the same telescope 
used by iDickman et al.l (jl992ft . iDickman et al.l mapped 
IC 443 with full beam spacing (i.e., undersampled), but 
we obtained fully sampled map using on-the-fTy (OTF) 
technique. The details of observations are described in 
§ O Our observations find multiple components of ambi- 
ent molecular clouds toward IC 443. The characteristics 
of these multiple cloud complexes are presented in § 03 In 
§ IH we try to identify ambient clouds that are physically 
associated with IC 443. The implication of our results 
on the environment of IC 443 is discussed in § [5] and § [6] 
summarizes the paper. 

2. OBSERVATIONS 
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A ~ 1° x 1° area around IC 443 was observed in 12 C0 
J=l-0 (115.2712 GHz) and HCO+ J=l-0 (89.18852 
GHz). All observations were obtained using the FCRAO 
14 m telescope during 2003 January. The single sideband 
32 element focal plane array receiver SEQUOIA was used 
in OTF mapping mode. The dual channel correlator 
(DCC) allowed us to observe two frequencies simulta- 
neously. An area of 60' x 62' centered at (06 h 17 m 00 s , 
22°34 / 00 // ) was mapped (Fig. [Q) with both channels of 
DCC set to 12 CO but at different ^lsr, one at ^lsr = 
—45 km s _1 and the other at 35 km s _1 . A smaller area of 
48'x61' centered at (06 h 17 m 26 s , 22°34'08") was mapped 
with two channels of DCC set to HCO+ and HCN at 
f lsr = — 15 kms -1 , respectively. In this paper, how- 
ever, we focus on the 12 CO and HCO + data as the inter- 
pretation of the HCN data is complicated due the com- 
plexity of its hyperfine structure. The FWHM beam- 
size of the FCRAO telescope is - 45" at 115 GHZ and 
- 55" at 89 GHz. The OTF data are regridded on a 
regular grid of 15" pixels. The DCC has a bandwidth 
of 50 MHz and 1024 spectral channels per IF channel, 
leading to velocity resolution after Hanning smoothing 
of 0.254 and 0.328 kms -1 at 115 GHz (CO) and 89 GHz 
(HCO + ), respectively. For 12 CO, the velocity coverages 
were —110 kms -1 to +20 kms -1 , and —30 kms -1 to 
+100 kms -1 for each channel of the DCC. Previous stud- 
ies detected emission extending up to +30 kms -1 in the 
positive veloci ty and down to — 100 kms -1 in the neg- 
ative velocity (jSnell et al.ll2005[ and references therein). 
Our observations may not provide adequate velocity cov- 
erage for this very high velocity gas. We find that the 
lack of baseline is severe in some locations for the posi- 
tive velocity DCC channel data, while negligible for the 
negative velocity DCC channel data. As the negative 
DCC channel data cover virtually all the interesting ve- 
locity range (see Section [3]), we primarily use the neg- 
ative channel data. For HCO + , the velocity coverage 
was —100 kms -1 to +70 kms -1 . The characteristic rms 
noise of each spectral channel is ~ 0.2 K for 12 CO, and 
~ 0.5 K for HCO + . All the plotted data in this paper are 
in antenna temperature (TjQ, which may be converted to 
brightness temperature by dividing by the main beam ef- 
ficiency which is 0.45 and 0.50 at 115 GHz and 90 GHz, 
respectively. In Fig. [TJ we overlay the areas we observed 
on top of the 21 cm radio continuum image of IC 443 
(LEE08). 



3. RESULTS 

The spectra of 12 CO and HCO + lines averaged over 
the entire observed field are shown in Fig. [2] Most of the 
CO emission arises from two narrow velocity components 
at approximately —3 kms -1 and +5 kms -1 , presumably 
from ambient molecular cloud. However very weak and 
broad emission from the shocked gas can also be seen. 
The average HCO + spectra are dominated by broad lines 
from the shocked gas. The average spectrum of 12 CO in 
Fig.[2]shows an artifact around ^lsr= —20 kms -1 , which 
is the velocity of atmospheric line and we attribute this 
artifact to an imperfect sky subtraction. This artifact 
is only seen in the average spectrum and has negligible 
effect on individual 12 CO spectra. 



3.1. Shocked Molecular Gas 

The shocked molec ular gas is chara cterized by its 
broad line wings (see iSnell et all I2005L and references 
therein). For ^lsr < — 10 kms -1 and ^lsr > +10 
kms - , the emission is confined spatially to small re- 
gions with linewidths that are very broad. This emission 
is assumed to arise almost entirely from the shocked gas. 
Due to projection effects, the velocity of some of the 
shocked gas is close to kms -1 (Dickma n et al.l fl992) 
and a significant fraction of the shocked gas is indeed in 
the velocity range between ^lsr= —10 to +10 kms -1 . 
But the study of this low velocity shocked gas with 12 CO 
is difficult due to the contamination by the ambient gas 
emission (Fig. [2j). HCO + emission, on the other hand, 
has less contamination from the ambient gas. The chan- 
nel maps and the position- velocity (PV) maps of the ob- 
served HCO + line are shown in Fig. [3] and in Fig. |4j In 
particular, the PV maps clearly show that most of the 
emission has a line width > 10 kms -1 , and thus is from 
the shocked gas. Some ambient gas is picked up as ab- 
sorption features and also as emission features. Thus 
HCO + is better suited for studying the overall distribu- 
tion of the shocked molecular gas. In Fig. HJ the distri- 
bution of high- velocity 12 CO emission is compared to the 
distribution of HCO + emission. We note that the 12 CO 
image at higher velocity can provide more sensitive and 
detailed distribution of shocked gas compared to HCO + 
data with a caveat that the lower velocity shocked gas is 
not included. 

The distribution of shocked molecular gas takes a 
shape of an incomplete ring, which consists of several 
bright clumps and fainter emission that seems to connect 
these clumps. The shocked gas emission is brighter in the 
south, while the emission in the north is weaker. There is 
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Figure 2. (a) Average 12 CO (J =1-0) spectra of entire field ob- 
served. Two spectra, one centered at ^lsr— — 45 kms -1 and the 
other at ^lsr— 35 kms -1 , are shown. The velocity range for chan- 
nel maps in Fig. [5] is marked. The bump around —25 kms -1 is 
an artifact due to an incomplete sky subtraction, (b) Same as (a) 
but y-range is adjusted to show the faint broad lines, (c) Average 
HCO+ spectra of entire field observed. Note that the temperature 
scale is different. 
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Figure 3. HCO + channel maps for velocities between —10 kms 1 and +10 kms 1 . The emission is predominantly from shocked molecular 
gas. However, there is some contribution from the ambient clouds, as can be identified by narrow lines in the PV map (Fig. 3}. 



a general tendency of the systematic velocity of shocked 
gas getting more positive with the increasing declination. 
This large-scale velocity gradient is consistent with pre- 
vious studies which suggested that the large-scale distri- 
bution of the shocked gas in IC 443 could be described 
as an expanding ring (jBurton et al.lfl988c iDickman et al.l 
Il992[ ). On the other hand, we note that long (> a few ar- 
cminutes) filaments of shocked gas are only prominent in 
the southeastern (SE ridge) and western parts (W ridge) 
of the proposed ring (Figure |4]) and the gas kinematics 
are dominated by rapid variation on the scale of < a few 
arcminutes, which are referred as shocked "clumps" . The 
PV plot in right ascension (the top panel of Fig. [4]) shows 
a tilt of the emission in clumps B and C. In clump B, the 
emission shifts to smaller R.A. at higher shock velocities, 
while for clump C, the emission shifts to larger R.A. at 
higher shock velocities. Little shift is seen in the PV 
plot in declination. The velocity sh i ft see n in Clump B 
was first detected bv IDickman et al.l (|1992l ) in both 12 CO 
and HCO + emission. They suggested that this shift was 
a sign that the gas was being ablated and accelerated 
away from the shocked clump. It is interesting that we 
see a similar effect in Clump C, although the ablated gas 
is moving off in the opposite direction. Although the ge- 
ometry of the interaction of the SNR with the ambient 
gas is complicated, these tilts in the PV diagram provide 
an indication of the direction of the shock. 

3.2. Ambient Molecular Gas 
The channel maps of 12 CO in the velocity range of 
^lsr= —10 kms -1 to +10 kms -1 , where most of the 
emission traces ambient clouds, are shown in Fig. [5l 
IC 443 is located toward the Galactic anticenter direc- 
tion, and no ambient molecular gas emission is found out- 
side of this velocity range. The channel maps in Fig. [5j 



together with the mean 12 CO spectrum in Fig. [2j sug- 
gest that there are two distinct complexes of molecu- 
lar gas, one in the negative velocity of —10 kms -1 < 
t?LSR < kms -1 and the other in the positive velocity of 
+3 kms -1 < ^lsr < +10 kms -1 . 

The clouds at negative velocities show significant 
channel-to-channel variation. To get a grasp of the over- 
all structure of the cloud complex, we generate an inte- 
grated temperature map, a peak temperature map and a 
first moment (intensity weighted velocity) map of 12 CO 
emission in the velocity range — 10 < ^lsr < kms -1 , 
which are shown in Fig. [6l Also, in Fig. [7J an RGB 
composite image is shown, where each red, green and 
blue channel corresponds to an integrated temperature 
map of different velocity ranges (see the figure caption 
for more details). The integrated temperature map pro- 
vides an overall view of the spatial distribution of this 
gas, although there is some contamination by emission 
from the shocked gas. Most of the ambient gas is dis- 
tributed along the northwest-southeast direction across 
the field, roughly corresponding to the interface between 
Shells A and B (Fig. [7j). The peak temperature map is 
similar to the integrated temperature map, however, sev- 
eral small clouds (~ 1') with relatively higher peak tem- 
peratures (> 5 K) stand out against the diffuse clouds 
whose peak temperatures are generally less than 3 K. 
These small clouds are even more prominent in the PV 
maps in Fig. [U Fig. [8] clearly reveals two distinguished 
kinematic characteristics of the small clouds compared 
to the diffuse clouds. First, the linewidths of these small 
clouds are much narrower than the diffuse clouds. Sec- 
ond, most of these small clouds show clear offset in ve- 
locity relative to the diffuse emission. Although some of 
these small clouds are located toward the diffuse clouds, 
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Figure 4. HCO + integrated temperature map on the lower-left panel, and position-velocity maps on the upper and right panels. The 
upper panel is an R. A. -velocity map of peak temperatures along the decl. axis. The right panel is a velocity-decl. map of peak temperatures 
along the R.A. axis. The values are derived from the whole velocity range. The HCO+ integrated temperature map has an rms error of 
0.3 Kkms -1 . The R.A. -velocity map and the velocity-decl. map has an rms error of 0.03 K. The contours represent 12 CO integrated 
temperature map over velocity ranges of —17 to —10 kms -1 and +1 to +3 kms -1 , chosen to minimize the contribution from the ambient 
gas. The contour levels are 1.2, 2.6, 6, 10 Kkms -1 and the rms error is 0.3 Kkms -1 . The shocked clumps are also indicated as in 
Fig.^b). Some HCO + emission is from th e am bient gas as can be identified from their narrow line widths. These ambient gas components 
are associated with the small clouds (see § 13. 3|) . 



a significant velocity offset between them suggests that 
most of the small clouds are not related to the diffuse 
clouds. Thus, these small clouds are considered to be 
separate clouds from the diffuse clouds. For the purpose 
of the discussion, we denote the diffuse clouds between 
—6 kms -1 and —1 kms -1 as "—3 kms -1 clouds" and 
the small high brightness temperature clouds as "SCs". 
The characteristics of the SCs will be further discussed 
in § 13.31 The first moment map mostly traces the ra- 
dial velocity variation within the —3 kms -1 clouds and 
suggests that the main body of the cloud has a small, 
systematic velocity shift from east to west. A similar 
tendency can be noticed in the RGB composite image 

(Fig. . 

In the positive velocity range, relatively fainter emis- 
sion from ambient clouds is seen between ^lsr = +3 and 
+10 kms -1 . We call these clouds as "+5 kms -1 clouds". 
Emission between kms -1 and +3 kms -1 is predomi- 
nantly from shocked molecular gas. Fig. [9] shows inte- 
grated temperature map of ^lsr = +4 to +10 kms -1 , 
which represent mostly the ambient gas with some con- 
tribution of the shocked gas. The clouds seem to be 
divided into two parts: the western part and the eastern 
part. The western part is brighter and shows a promi- 
nent filamentary structure. The eastern part is fainter 



and shows a number of small clouds. It is noteworthy 
that the western and the eastern parts are separated by 
Shell A (Fig. [9j). 

3.3. Small Clouds (SCs) 

The 12 CO peak temperature map in Fig. [6] and the 
PV maps in Fig. reveal the existence of several SCs 
(~ 1 ') with relatively high brightness temperature and 
narrow line width. To identify these clouds in more 
systematic way, we utilize two characteristic quantities: 
the peak temperatures (T* e3k ) and the equivalent line 
width (A^eq). Av eq is defined as the ratio of the in- 
tegrated temperatures 7^* t to Tp eak . These tempera- 
tures are measured in the velocity range of —10 kms -1 
to kms -1 . Fig. [10] shows a two dimensional his- 
togram of these two quantities measured from the pix- 
els in the observed area. Most pixels have T* eak < 4 K 

and 1 < v eq < 4 kms -1 , and we find that they are 
mostly associated with —3 kms - clouds. On the other 
hand, there are pixels with relatively high peak temper- 
ature (T* eak > 4.5 K) with relatively narrower equivalent 

line width (v eq < 3 kms -1 ). We find these pixels cor- 
respond to the SCs as suspected. Based on the two- 
dimensional histogram, we adopt a criteria of T* eak > 
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Figure 5. 12 CO channel maps for velocities between —10 kms 1 and +10 kms 1 , mostly tracing the ambient molecular gas. Each panel 
represents mean temperature for Av = 0.5 kms -1 , where their central velocity is shown in the top-left corner. 
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Figure 6. (a) 12 CO integrated temperature map in velocity range 
of —10 to kms -1 , mostly tracing the ambient gas. (b) 12 CO peak 
temperature map over the same velocity interval. The numbers 
denote the location of small clouds (SCs) discussed in § 3.3. (c) 
The intensity weighted mean velocity (first moment) of the 12 CO 
emission. The contour s ma rk the regions which satisfy our selection 
criteria for SCs (see § I3.3l for details). 



4.6 K and v eq < 3.5 kms -1 , which results in the identi- 
fication of 12 SCs, and these are identified in Fig. [6fb) 
and labeled with numbers (e.g., SC 01). The selection 
process includes all of the SCs that would have been cho- 
sen by eye, as we are primarily selecting them based on 
peak temperature. 

The 12 CO line profiles of the 12 SCs are shown in 
Fig. [ITJ These are average spectra extracted from the 
contours shown in Fig. Fig. Efb), and they all show nar- 
row lines with high peak temperatures as expected. We 
fit the spectra with Gaussian profiles and the results are 
summarized in Table [TJ Often, the spectra need mul- 
tiple Gaussian components to be adequately fit. The 
additional components are contributed by emission from 
either the diffuse clouds or the shocked gas. 



Central velocities of these SCs range between —3 and 
—8 kms -1 , with a median of ~ —6 kms -1 . Other than 
SC 01, 07/08 (see the note in Table 1), they have peak 
velocities more negative than —5 kms -1 , and are well 
separated from —3 kms -1 clouds in the velocity space 
(Fig. EJ. Their line widths are between 1 and 2 kms -1 . 
The HCO + profiles are also plotted in Fig.fTTIwhen avail- 
able. The SCs tend to have similar narrow emission lines 
in HCO + . The 12 CO maps of ~ 8' x 8' regions around 
individual SCs at their peak intensity channels are shown 
in Fig. [121 Part of the emission seen around SCs 05, 06, 
07/08, 10, and 11 is due to shocked gas, and the associ- 
ation of these SCs to the shocked molecular gas will be 
further discussed in § 14.11 

4. IDENTIFICATION OF AMBIENT CLOUDS 
INTERACTING WITH IC 443 

The shocked molecular gas in IC 443 is mostly found 
within Shell A. No s hocked molecular gas has been identi- 
fied outside Shell A (jDickman et al.lfl992h . even with our 
new wide-field, fully-sampled mapping observations. The 
signature of the shocked molecular gas is pr ominent along 
the s outhern boundary of Shell A (e.g., Bu rton et al.l 
Il988f ). On the other hand, the nature of the ambi- 
ent molecular clouds that are associated with shocked 
molecular gas in IC 443 has not been clearly understood. 
Our observations reveal the existence of three different 
cloud complexes (—3 kms -1 clouds, SCs and +5 kms -1 
clouds) toward IC 443. Here we investigate possible as- 
sociation of these cloud complexes with IC 443. 

4.1. SCs 

The SCs around IC 443 are identified for the first time 
in our observations. We find that their spatial distribu- 
tion is closely related to the distribution of shocked gas in 
IC 443. In Fig.[l3j we overlay the locations of the SCs on 
top of the distribution of the shocked molecular gas. The 
positions of SCs 05, 06, 10, and 11 are coincident with 
those of the shocked molecular clumps. Also, the south- 
eastern tip of SC 03 is coincident with Clump G, and the 
southwestern tip of SCs 07/08 with Clump E (we note 
that these SCs are spatially more extended than those 
regions defined by our selection criteria). While some of 
the positional coincidence could be by chance, the fact 
that a large fraction of the SCs shows their projected po- 
sitions on the sky being coincident with those of shocked 
clumps suggests that they are more likely physically as- 
sociated. 

Evidence of the physical association of the SCs with 
the shocked molecular gas is best seen in SC 05. Fig. [T4l 
shows PV maps of 12 CO and HCO+ along an NE-SW 
direction and the location of this PV cut is indicated in 
Fig. [T3l The spatial extent of the shocked gas is nearly 
identical to that of SC 05, unlike the emission from the 
more spatially extensive —3 kms -1 clouds. This coinci- 
dence strongly suggests that SC 05 is physically related 
to the shocked gas. 

We further find that the locations of the three OH 1720 
MHz maser sources currently known in IC 443 are spa- 
tially close to the SCs or even coincident (Fig. [12]). Fur- 
thermore, the systematic velocity of the maser source in 
SC 05 (-6.1 ±0.1 kms -1 . lHewitt et al.ll2006l ) is coinci- 
dent with the systematic velocity of SC 05 itself. The ve- 



8 



Lee, J. -J. et al. 




6 h 19 m 6 h 18 m 6 h 17 m 6 h 16 m 6 h 15 m 

<^2000 

Figure 7. RGB composite image where each R, G, and B channel corresponds to 12 CO integrated temperature maps of velocity ranges 
—0.5 to —3.4, —3.4 to —5.1 and —5.1 to —7.0 kms -1 , respectively. The contours are 21cm radio continuum image from LEE08. The 
contour levels are 10, 20, 40, 60, 80 mJy. 



Table 1 

Spectral parameters of the small clouds (SCs; see § 13. 3|) 



Component 1 Component 2 



SC No. 


a (2000) 


S (2000) 


T* [Kj 


v 






%WHM 


T* [Kj 


v 






^FWHM 


01 


06 h 15 m 20?l 


22°53 / 43 // 


5.7 ± 0.5 


-3.1 ± 


0. 


1 


1.6 ± 0.1 


0.9 ± 0.3 


-4.6 ± 


0. 


5 


2.7 ± 0.6 


02 


06 h 16 m 14!4 


22°45 / 14 // 


6.0 ± 0.1 


-7.5 ± 


0. 


1 


0.9 ± 0.1 












03 


06 h 16 m 29!6 


22°36 / 14 // 


7.3 ± 0.1 


-6.2 ± 


0. 


1 


1.2 ± 0.1 


2.6 db 0.1 


-4.1 ± 


0. 


1 


3.6 ± 0.1 


04 


06 h 16 m 25!3 


22°22 / 44 // 


6.6 ± 0.1 


-6.0 ± 


0. 


1 


1.2 ± 0.1 


0.5 ± 0.1 


-3.5 ± 


0. 


1 


2.0 ± 0.2 


05 


06 h 17 m 29!l 


22°22 / 44 // 


6.7 ± 0.1 


-6.2 ± 


0. 


1 


1.3 ± 0.1 












06 


06 h 18 m 08!l 


22°27 / 44 // 


4.4 ± 0.2 


-7.9 ± 


0. 


1 


1.5 ± 0.1 


0.9 ± 0.1 


-9.2 ± 


0. 


2 


5.2 ± 0.3 


07 


06 h 18 m 08 s 2 


22°34 / 58 // 


4.1 ± 0.3 


-4.1 ± 


0. 


1 


2.2 ± 0.1 


1.6 ± 0.1 


-4.7 ± 


0. 


1 


5.4 ± 0.3 


08 


06 h 18 m 17!9 


22°36 / 28 // 


6.6 ± 0.1 


-4.6 ± 


0. 


1 


2.4 db 0.1 












09 


06 h 19 m 04!6 


22°43 / 57 // 


3.5 ± 0.4 


-6.7 ± 


0. 


1 


0.7 ± 0.1 


2.1 ± 0.2 


-7.1 ± 


0. 


1 


1.7 db 0.1 


10 


06 h 17 m 48 s 7 


22°40 / 14 // 


4.6 ± 0.2 


-5.3 ± 


0. 


1 


0.7 ± 0.1 


1.0 ± 0.1 


-6.4 ± 


0. 


1 


7.0 ± 0.3 


11 


06 h 17 m 48?7 


22°41 / 59 // 


4.0 ± 0.5 


-6.0 ± 


0. 


1 


0.9 ± 0.1 


1.1 ± 0.2 


-6.4 ± 


0. 


1 


2.3 ± 0.3 


12 


06 h 17 m 06!5 


23°01 / 45 // 


5.0 ± 0.1 


-6.0 ± 


0. 


1 


1.4 ± 0.1 


1.1 ± 0.1 


-3.7 ± 


0. 


1 


2.1 ± 0.1 



Note. — SCs 07 and 08 are positionally close and have similar systematic velocities, and they will be labeled 
together as SCs 07/08 when necessary. 



locity of the maser source near SC 06 (—6.9 ±0.1 kms - , 
iHewitt et al.ll2006[ ) is also in agreement with the velocity 
of SC 06. The velocity of the source near SC 03 (the 
source that is coin cident with Clump G) is —4.55 kms -1 
([Hewitt et al.ll2QQ6D . While this is slightly different from 
the peak velocity of SC 03, it is still within the velocity 
wing of SC 03. We consider that this maser source could 
be associate d with SC 03. 

Recently, llndriolo et al.l ((2010) searched for H3" ab- 
sorption features in six sight lines towards IC 443. From 



the absorption features seen in two sight lines, high ion- 
ization rates are inferred which are attributed to the in- 
creased cosmic ray ionization rate near IC 443. The two 
sight lines with high ionization rates are indeed located 
close to the SCs (see Fig. [12]). The centroid velocities of 
the absorption features are between —8 and —6 kms -1 , 
being consistent with the velocity range of the SCs. 

Therefore, we found that a good fraction of the SCs 
show direct and/or indirect association with the shocked 
molecular clumps. We propose that the SCs are at the 
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Figure 8. 12 CO peak temperatures map (identical to Fig.[6fb)) together with position-velocity maps. The upper panel is a R. A. -velocity 
map of peak temperatures along the decl. axis. The right panel is a velocity-decl. map of peak temperatures along the R.A. axis. The 
images are shown in a squared scale to better visualize the locations of the SCs. The numbers are identification numbers for SCs. The 
central velocities of SCs are adopted from vq of Component 1 in Table [T] 
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Figure 9. 12 CO distribution of +5 kms 1 clouds, integrated over 
velocity range of +4 kms -1 to +10 kms -1 . The contours are 12 CO 
integrated temperature map as in Figure[4] tracing the distribution 
of shocked molecules. 

same distance as IC 443, and that some of them are in 
direct interaction with IC 443. 

4.2. +5 kms -1 clouds 

The detailed morphology of +5 kms -1 clouds is also 
newly revealed in our observations. The overall morphol- 
ogy of the clouds, which is shown in Fig. [9j is intriguing. 




10 12 



Figure 10. (top) distribution of T* eak vs. Ai; eq (equivalent line 
width), where A^eq, defined as T* nt / T* 



have T peak~ 4 K and Av eq ~ 

T p * eak > 5 K and A^ eq < 3 kms" 



int/ -peak" While mOSt P ixels 

2 kms -1 , there are pixels with 
1 which are associated with SCs. 



The gray lines mark the selection criteria used to identify SCs. 
(bottom) Histogram of T* eak . A clear bump around 5 K is seen. 

The gray vertical line marks T* eak =4.6 K above which are neces- 
sary condition to be identified as SCs. 
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Figure 11. 12 CO line profiles of the SCs. The crosses indicate the line center and the full-width-at-half-maximum when fit by multiple 
Gaussian components (see Table [T}. HCO + spectra at the same locations are also shown as gray lines. No HCO + spectrum is available 
for SC 01 as the observed area for HCO + is smaller than that for 12 CO. 



As has been previously described, the clouds can be di- 
vided into two parts. The brighter western part has a 
prominent filamentary morphology. It extends from west 
to east. Its eastern boundary is sharp and we see a fila- 
ment of shocked gas along this boundary (clump G and 
W ridge). In Fig. [15j we compare the morphology of 
+5 kms" clouds to that of 90 fim far infrared emis- 
sion from IC 443. The far infrared emission, taken by 
AKARI telescope, is presumably from the shock heated 
dust and defines the outer boundary of Shell A. Interest- 
ingly, the +5 kms -1 clouds (both the western and the 
eastern parts) are only visible along and/or outside Shell 
A. Therefore, there is a good morphological coincidence 
between the +5 kms -1 clouds and the shock tracers. 

On the other hand, the systematic velocity of the 
shocked gas in Clump G and the W ridge is around 
—5 kms -1 , and thus is considerable offset from that 
of the +5 kms -1 clouds. The symmetric line profile 
of Clump G suggests that the shock direction is mostly 
perpendicular to our line-of-sight, thus shock geometry 



is unlikely to explain the 10 kms -1 difference in radial 
velocity. Can these kinematically distinct clouds be ad- 
jacent in space? IC 443 is likely a co re-collapse SNR 
located in the Gem OBI association ([Braun fc Stroml 
1986), and there could be clouds with different velocities 
due to the complicated kinematics within this OB asso- 
ciation. However, the difference of 10 kms -1 seems to 
require rather extraordinary condition. Therefore, while 
the morphology of the +5 kms -1 clouds is rather sugges- 
tive of its association with IC 443, there is no kinematic 
evidence that these +5 kms -1 clouds are responsible for 
any of the shocked molecular gas. We do not consider 
that the physical association of the +5 kms -1 clouds 
with IC 443 is conclusive. 

4.3. —3 kms -1 clouds 

The —3 kms -1 clouds show elongated morphology 
in the NW-SE direction, along the interface between 
Shells A and B. There is supporting evidence that 
most of thes e clouds are locat ed on the front side of 
the remnant (iTroja et al.ll2006[ and references therein), 
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Figure 12. 8 ; x 8' close-up view of each SC at its peak velocity (see Table Each panel is labeled with the SC number and the channel 
velocity used. The white contours are T* = 5 K, 7.5 K and 10 K. The plus signs near SCs 03, 05 and 06 are positions of OH maser sources 
and the star signs near SC 02 and 06 are sight lines along which absorption features are detected (see § 14.10 . The locations of SC 07 
and SC 08 are also marked with arrows for clarification. 



and the self-absorption f eatures in 12 CO and HCO + 
(|van Dishoec k et al.l IT993) are consistent with this con- 
clusion. Thus, any interaction of IC 443 with these clouds 
would likely produce blueshifted gas; however, we see 
both redshifted and blueshifted shocked gas. Also, we 
see little correlation between the locations of the shocked 
gas and the spatial distribution of the —3 kms -1 clouds. 
While the velocity gradient within the —3 kms -1 clouds 
might be due to some external perturbation, we could 
not find any hint that this is due to the interaction with 
IC 443 or its progenitor star. An OH 1720 MHz maser 
source associated with Clump G has a systematic velocity 
of —4.55 kms - (He witt et al.ll2006h . comparable to the 
velocity of —3 kms -1 clouds. However, the location and 
the velocity of the maser source are also consistent with 
those of SC 03, so the association of this maser source 



with the —3 kms clouds is considered ambiguous. We 
find no strong evidence for a spatial or kinematic connec- 
tion between the shocked gas and the —3 kms -1 clouds, 
thus we conclude that the remnant is not currently in- 
teracting with the —3 kms -1 clouds. 

5. DISCUSSION 

5.1. Origin of the Small Clouds & Environment of 
IC 443 

The physical characteristics of the SCs, size of ~ 0.5 
pc (~ V at 1.5 kpc) and velocity dispersion of ~ 
1 kms -1 , resemble sma ll clumps in molecular clouds (cf., 
iBergin fc Tafa lla 2007). We estimate the hydrogen col- 
umn density of SCs using the canonical conversion factor 
x (= N U J fT( 12 CO)dv) of - 2 x 10 20 cm- 2 K- 1 km- 1 s 



12 



Lee, J. -J. et al. 









02 -^s. 


22° 45' 


11 

% 

g 


10 


O 


§ 




* 1 

3 




22° 30' 










06 -»* 

N 


m 

05 


f 

Q «-04 


22° 15' 










6 h lg m 


6 h 17 m 

^2000 





Figure 13. Locations of the SCs overlaid over the average tem- 
perature map of HCO + . The gray rectangle is the extraction area 
for the PV map in Fig. [TH 
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Distance ['] 

Figure 14. PV map of 12 CO emission around SC 05 extracted 
from the slit in Fig. 1131 The distance is measured from the location 
of SC 05 and the value increases with increasing right ascension. 
The contours represent a PV map of HCO+ extracted from the 
same slit. The location of SC 05 in the PV map (distance ~ ; , 
i>lsr ~ — 6 km s ~ 1 ) is marked. 

fe.g. JDame et al.ll2QQlf ). and list these in Table [2 If we 
assume a path length of ~ 0.5 pc (1', which is a typical 
angular diameter of the SCs, at a distance of 1.5 kpc), we 
obtain average molecular hydrogen densities of the order 
of 10 3 cm -3 (higher densities must be present to produce 
the HCO + emission seen in most of the SCs, suggesting 
that the clouds could be clumped and have a small filling 
factor within the beam). The masses of the individual 
clouds are estimated to be of order of 10 M (Table [2]). 
This is comparabl e to the mass of in dividual shocked 
clumps (5-40 iDickman et al.l I1992T) . and also com- 
parab le to that of small clumps (cf., iBergin fc Tafalla 
l2007h . 

One of the key characteristics of the SCs is that they 
lack surrounding clouds. We propose that the SCs are 
analogous to the cometary globules, i.e., remnant cores of 
molecular clouds that have been exposed as their lower 
density gas envelopes have been removed by the prior 
activity of the SN progenitor star or nearby OB stars 
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Figure 15. Far IR 90 /im image taken with the AKARI satellite 
shown in gray scale. The green contours show the distribution of 
+5 kms -1 clouds (the gray scale in Fig. [9}. The blue contours 
show locations of SCs. The solid and dashed circles represent the 
location of 7-ray sources detected by MAGIC and VERITAS, re- 
spectively. 



Table 2 

H2 column densities and masses of SCs. 



SC# H 2 


Column Density 


Mass 




[xlO 21 cm" 2 ] 


[Mq] 


01 


4.3 ±0.5 


12.6 ± 1.6 


02 


2.5 ±0.1 


3.5 ±0.1 


03 


4.0 ±0.1 


57.7 ±0.9 


04 


3.6 ±0.1 


7.0 ±0.1 


05 


4.1 ±0.1 


13.2 ±0.2 


06 


3.0 ±0.2 


2.0 ±0.1 


07 


4.2 ±0.4 


2.7 ±0.3 


08 


7.4 ±0.1 


85.8 ±0.7 


09 


1.2 ±0.2 


1.2 ±0.2 


10 


1.5 ±0.1 


0.8 ±0.1 


11 


1.7 ±0.3 


0.4 ±0.1 


12 


3.2 ±0.5 


36.9 ±0.7 


Note. 


— H2 column densities 



are estimated from the fit parameters 
of component 1 in Table [T] We cor- 
rected for the main beam efficiency and 
adopted the canonical conversion factor 
X (= 7Vh 2 //T( 12 CO)^) of ~ 2 x 
10 20 cm- 2 K- 1 km- 1 s (e.g., IDame et all 
[20011 V To estimate the mass, we multiply 
the H2 mass by a factor of 1.36 to account 
for the mass of helium. The errors in this 
table are formal errors from the spectral fit 
and do not include any systematic error. 

(jReipurthl fl983h . For SCs 01 and 12, located farther 
from IC 443, the removal of lower density envelopes may 
have not been complete as we see some diffuse ambient 
gas around them. Some of the SCs are now being im- 
pacted by the SNR and the gas is being ablated and ac- 
celerated producing the shocked molecular clumps (e.g., 
SC 05 and the associated shocked gas). The interaction 
of small clouds with shock waves have been comprehen- 
sively studied by numerical simulations (e.g. JKlein et al.l 
1994). Upon interaction, transmitted shocks will propa- 
gate into the cloud core and the envelope will be stripped 
away. This nicely explains the spatio-kinematic structure 



Identification of Ambient Molecular Clouds around IC 443 



13 



22° 30' 



Gray : 12C0 Tpeak(-9.2 5.4 km/s) 



22° 25' 




22° 20' 



22° 15 



Tmean(-10 
^> Tmean(-15 
Cl"> Tmean(-35 



15 km/s) 
35 km/s) 
70 km/s) 



gh 18 m QQS 



40 s 



20 s 



gh 17 m QQS 



<^2000 



Figure 16. Gray scale shows a peak temperature map of 12 CO in velocity range —0.2 and —5.4 kms -1 , mostly representing the ambient 
gas. Red, green, and blue contours are integrated temperatures of 12 CO for velocity range of —10 to —15 kms -1 , —15 to —35 kms -1 , and 
—35 to —70 kms -1 , respectively. They show emission from shocked molecular gas. The long rectangles are extraction regions for the PV 
maps in Fig. IT71 
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Figure 17. PV maps along the slit shown in Fig. 1161 The upper 
panel is for 12 CO, and the lower panel is for HCO + . The gray 
contours, both in upper and lower panels, are those of HCO + . The 
offsets are distance in arcminute measured from Clump C (marked 
as a plus sign in Fig. I16|) and the offset increases with increasing 
right ascension. 



of the gas around SC 05. In this picture, some part of the 
core is already shocked and ablated, while the remaining 
part is yet to be shocked. The ambient gas of SC 05 itself 
is unperturbed part of the core, and the high velocity gas 
would represent a flow of gas that is primarily ablated 
from the core. 

In figure [171 we show a PV diagram along a slit extend- 
ing NW-SE through SC 05 and covers Clumps B and C 
and SC 05 (see Figure [16] for the location of the slit). 
It indicates that the shocked gas associated with SC 05 
is quite similar to that of Clumps B and C. These two 
clumps show more negative velocity emission than that 
of SC 05, but no SC-like ambient cloud is detected for 
them. Thus, we suggest that Clump B and C represent 
the later evolutionary stage of the shock-core interaction 
and the cores associated with Clumps B and C are al- 
ready been destroyed. On the other hand, some of the 
SCs (e.g., SC 04) do not have associated shocked gas. A 
compact ambient cloud similar to the SCs is also seen 
near the slit offset of ~ 2 1 in Figure [TTJ Its characteris- 
tics are similar to that of other SCs, except that the peak 
temperature is not high enough to meet our SC criteria. 
We suggest that this cloud is yet to be impacted by the 
SNR shock. 

The absence of small SC-like features with Clumps B 
and C points toward the rapid destruction of the core 
following the passage of the shock. The clouds are likely 
accelerated and destroyed within several cloud crushing 
time intervals, where the cloud crushing time is defined as 
t C c = X 1 a o/ v b and ao, each represents the ratio of 
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the density of the cloud to that of the intercloud medium, 
the cloud size, and the sho ck velocity in the intercloud 
medium (Kle in et al J 11994 ). If we adopt ao of 1 pc and 
X ~ 100, we obtain r cc ~ 1 x 10 5 (^/lOOkms -1 ) -1 
yrs. Vb is uncertain, but recent X-ray observations 
(jTroja et al.l l2006) seem to suggest ~ 500 kms -1 while 
other observations tend to give lower velocities (e.g., 
see lChevalierlll999f ). Adopting ~ 500 kms -1 yields 
r cc ~ 2 x 10 4 yrs, which is c omparable to th e estimated 
remnant age (~ 30,000 vrs. lChevalier|[l999[ ). Hence the 
shock crushing may not be an efficient mechanism to 
destroy SCs. However, the i ntercloud shock in IC 443 
could be in radiative phase (Ch evalierl [1999) and the 
large postshock density of the intercloud shock may have 
significantly accelerated the destruction of the clouds. 
The SCs are also subject to evaporation. The evapora- 
tion time for classical thermal conductivity is given by 

W = 3.3 x 10 20 n c (^) 2 77 5/2 (ICowie fc McKedfl977h . 
where n c is the hydrogen number density of the cloud 
and Ti is the temperature of the shocked int ercloud gas. 
Adop ting n c ~ 10 3 cm -3 and % = 4 x 10 6 K (jTroia et all 
2006) yields t evap > 2 x 10 6 yrs. This is much greater 
than the remnant age so the evaporation process should 
not have a significant effect on the evolution o f the SCs. 

On the other hand, as has been discussed in I Chevalierl 
(1999), O-type stars likely clear out molecular clouds far- 
ther than a 15 parsec in radius. Thus the existence of 
the SCs suggests that the progenitor star of IC 443 would 
not have been an O-type star. It would be more likely an 
early B star (initial mass of 8-12 M ) whose photoion- 
izing radiation and wind power are enough to clear out 
diffuse gas but not enough to clear out dense cloud cores. 

5.2. j -ray Emission 

IC 443 is one of the SNRs known to be associated 
with 7-ray sources . Recent observ ations with MAGIC 
(| Albert et al.ll2007l ) and VERITAS (jAcciari et al.ll2009l ) 
revealed a TeV source in IC 443. The emission is seen 
near Clump G (Fig. [T5]) . The origin of the 7-ray emis- 
sion in SNRs is currently not fully understood. While 
it is likely from the hadronic interaction of cosmic ray 
protons with dense ambient medium. One of the keys 
to understanding the origin of the 7-ray emission is the 
characteristics of the ambient clouds that serve as target 
particles of hadronic interactions. 

The location of the detected TeV 7-ray emission is close 
to the location of the shocked clump G and this is a region 
where all three components of the ambient molecular 
clouds are seen. It is not clear which cloud(s) is the target 
of cosmic ray protons. The fact that the 7-ray emission is 
only detected in this region implies that the target parti- 
cles are particularly abundant in this area. In this regard, 
—3 km s _1 clouds could not be the primary target clouds. 
The recent VERITAS observation suggests that the 7- 
ray emission is elongated in the north-south direction 
(jAcciari et al.ll200H ) and +5 kms -1 clouds cloud could 
provide major target particles. The intensity and spec- 
tral shape of the 7-ray emission depend on the distance 
and m ass of the target clouds. Recently, iTorres et al.l 
(|2010f ) found that the observed 7-ray spectra of IC 443 
could be explained by two cloud components; a cloud of 
~ 350 Mq close to the remnant (~ 4 pc) and a more 



distant (~ 10 pc) massive one of ~ 4000 M . The closer 
cloud component may correspond to the SC 03 or the 
shocked clump G while the distant one would be the 
+5 kms -1 clouds or the —3 kms -1 clouds. 

6. SUMMARY 

In this paper, we present fully sampled, spectroscopic 
imaging of 12 CO and HCO + molecular line emission to- 
ward the supernova remnant IC 443 to investigate its 
molecular environment and to survey the area for ad- 
ditional zones of shock emission. While all previously 
known shocked clumps are detected, no additional re- 
gions of shocks are identified in our observations. The 
ambient molecular gas toward IC 443 shows three kine- 
matically and morphologically distinct components. The 
—3 kms -1 clouds are elongated in the NW-SE direction 
along the interface between Shells A and B of IC 443. 
However, there are no signposts of interaction between 
IC 443 and this molecular component. A fainter cloud 
at +6 kms -1 is distributed along a northeast-southwest 
axis. While the morphology of its western edge is sugges- 
tive of an interaction, the gas velocities are too displaced 
from those of the shocked clumps to firmly establish its 
connection to the remnant. Twelve small clouds (SCs) 
with narrow line widths and high peak temperature are 
spatially coincident with shocked molecular clumps and 
other tracers of interaction with the remnant. We pro- 
posed that these SCs are the residual cores of the parent 
molecular cloud that have been exposed due to the radi- 
ation field and stellar winds of the progenitor of IC 443 
and nearby OB stars. IC 443 is currently encountering 
some of these remaining cores as indicated by the asso- 
ciation of the SCs with shocked molecular clumps. Some 
shocked molecular clumps are not associated with SCs. 
In these cases, the cores have been recently destroyed 
by the shock interaction. Finally, we note that all three 
cloud components are found toward the sight line of the 
7-ray source. These components likely provide the reser- 
voir of protons needed to account for this 7-ray source. 

We thank to the anonymous referee for valuable com- 
ments. This work is supported by grant AST-0838222 
from the National Science Foundation. 
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